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Abstract. We study the dynamics of Rydberg ions trapped in a linear Paul trap, 
and discuss the properties of ionic Rydberg states in the presence of the static and 
time-dependent electric fields constituting the trap. The interactions in a system of 
many ions are investigated and coupled equations of the internal electronic states and 
the external oscillator modes of a linear ion chain are derived. We show that strong 
dipole-dipole interactions among the ions can be achieved by microwave dressing fields. 
Using low-angular momentum states with large quantum defect the internal dynamics 
can be mapped onto an effective spin model of a pair of dressed Rydberg states that 
describes the dynamics of Rydberg excitations in the ion crystal. We demonstrate that 
excitation transfer through the ion chain can be achieved on a nanosecond timescale 
and discuss the implementation of a fast two-qubit gate in the ion chain. 
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Figure 1. Typical length scales in a chain of cold Rydberg ions in a linear Paul trap. 
The external trapping frequency is in the order of MHz with a corresponding oscillator 
length Xho of approximately 10 nm. The intcrparticlc spacing £, set by the equilibrium 
between the Coulomb forces among the ions and the external confinement, is typically 
about 5 /im. The third length scale is the size of the Rydberg orbit OR. y . Due to the 
scaling proportional to the square of the principal quantum number n it can assume 
values in the order of 100 nm and therefore become significantly larger than echo- In 
this regime the Rydberg ion cannot be considered as a point particle but rather as a 
composite object, and its internal structure must be taken into account. 



1. Introduction 

Rydberg states correspond to the infinite series of excited bound states in a Coulomb 
potential with large principal quantum number n. In view of their "macroscopic" size, 
a Ry ~ n 2 a Q with oq the atomic Bohr radius, Rydberg states have remarkable properties, 
as reflected, for example, in their response to external static and time dependent electric 
and magnetic fields [lj. While the single particle physics of Rydberg atoms has been 
the subject of intensive studies in the context of laser spectroscopy, recent interest 
has focused on exploiting the large and long-range interactions between laser excited 
Rydberg atoms to manipulate the many-body properties of cold atomic ensembles. 
Examples include recent seminal experiments on frozen Rydberg gases obtained by 
laser excitation from cold atomic gases, demonstrating in particular a dipole-blockade 
mechanism [21 El HI El El II] , which in sufficiently dense gases prevents the excitation of 
ground state atoms in the vicinity of a Rydberg atom, and proposals for fast two-qubit 
quantum gates between pairs of atoms in optical lattices [8]. Furthermore, Rydberg 
atoms have been proposed to serve as model systems for studying coherent transport 
of excitations [9] - a mechanism which is of great importance for coherent energy 
transfer in biological systems, e.g. in light-harvesting complexes [TOl [TT]. While these 
investigations have so far concentrated on neutral atoms, we are interested below in 
describing the properties of laser excited Rydberg ions stored in a Paul trap, in particular 
the interplay between trapping fields and Rydberg excitations, and the associated many- 
body interactions in a chain of cold trapped Rydberg ions. 

Atomic ions can be stored in static and radio frequency (RF) electric quadrupole 
fields constituting a Paul trap [12], where for sufficiently low temperature they form 
a Wigner crystal [131 [Til, [T5l [16] . Using laser cooling the ions can be prepared in the 
vibrational ground states of the phonon modes of the crystal. Internal electronic states 
of the ions can be manipulated with laser light and entangled via the collective phonon 
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modes. Here we consider a situation where ions initially prepared in their electronic 
ground state are excited with laser light to high lying Rydberg states. In contrast to 
ions in low lying electronic states, a Rydberg ion in a Paul trap must be understood as 
a composite object, where the Rydberg electron is bound by the Coulomb force to the 
doubly charged ion core, but both the Rydberg electron and the ion core move in the 
electric fields constituting the trapping potentials. This is reminiscent to the situation 
in which Rydberg atoms are confined by very tight magnetic traps [T8l [T9] . We are 
particularly interested in the parameter regime where the size of the Rydberg orbit 
<2R y is larger than the localization length of the doubly charged ion core x^ around its 
equilibrium position, but still much smaller than the average distance ( between the 
ions in the Wigner crystal, i.e. Xho < ^Ry <^ C ( see Fig- CO)- O ur goal is to provide a 
description of the interaction and quantum dynamics of such a ID string of Rydberg 
ions in a linear ion trap. 

In comparison to neutral Rydberg atoms, a number of features and differences 
emerges for Rydberg ions in a Paul trap, which will be illuminated in the present work. 
First of all, the character of a Rydberg ion as a composite object gives rise to an intrinsic 
coupling of electronic and external motion in the presence of the electric trapping fields. 
This will be shown to result in renormalized trapping frequencies for Rydberg ions 
compared to their ground state counterparts. Furthermore, the interaction among ions 
is not - as in the neutral case - governed by the dipole-dipole force alone but also by 
the charge- charge, dipole-charge and charge-quadrupole interaction. The interplay of 
this rich variety of interactions among the ions and the external trapping fields will be 
analyzed. 

In our study we focus on ionic Rydberg states with low angular momentum quantum 
number and correspondingly large quantum defect, and large fine structure splitting, 
as these states can be most simply described and most easily excited by laser light 
from electronic ground states. For typical ion trap parameters [12J static dipolar and 
van-der-Waals interaction among ions in these Rydberg states is shown to be small 
compared to the energy scale set by the external trapping frequencies of the ions. 
Thus, in order to establish substantially stronger interactions, we employ additional 
microwave (MW) fields driving transitions between Rydberg states. This leads to large 
oscillating dipole moments, which result in remarkably strong controllable dipole-dipole 
interactions between the ions. Our findings show that for strong interactions the internal 
and external dynamics of the ion chain approximately decouple such that a "frozen" 
Rydberg gas is formed. In this limit the Hamiltonian describing the electronic dynamics 
can be formulated in an effective spin- 1/2 representation and involves resonant and off- 
resonant dipole-dipole interaction terms with coupling strengths of the order of several 
hundred MHz. Based on this effective spin model we demonstrate resonant excitation 
transfer on a nanosecond timescale from one end to the other of the ion chain. Moreover, 
we show that a two-qubit conditional phase gate between adjacent ions, based on the 
dipole-dipole interaction, can be realized on a time scale, which is much shorter than 
both the external dynamics and the radiative lifetime of the involved Rydberg states. 
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The paper is organized as follows: In Sec.[2]we derive the Hamiltonian of the coupled 
internal and external dynamics of a single Rydberg ion in a linear Paul trap. We analyze 
the effects of the trapping fields onto the electronic properties of an ion excited to a 
Rydberg state, introduce MW-dressed Rydberg states and study the coupled electronic 
and external equations of motion within a Born-Oppenheimer approach. In Sec. [3] we 
turn to the analysis of several trapped Rydberg ions, discuss the variety of interactions 
among the ions and derive the effective spin- 1/2 Hamiltonian for the electronic dynamics 
in the ion chain. Finally we illustrate the quantum dynamics of the system by studying 
resonant excitation transfer and present the fast two-qubit gate scheme. 

2. Single Trapped Rydberg Ions as a Composite Object 

2.1. Single Rydberg Ion in a Linear Paul Trap 

A combination of static and time- dependent electric fields can be employed to confine 
charged particles in a restricted region of space. The electric potential of a quadrupole 
field of a standard Paul trap can be written in the form 



where u is the radio frequency (RF) drive frequency. The electric field gradients 
a and (3 are determined by the actual geometry of the experimental setup. In 
present experiments, typical parameters are a ~ 10 9 V/m 2 , (5 ~ 10 7 V/m 2 and 
uj = 2tt x 10.. .100 MHz (for details see e.g. Ref. [12]). 

We are interested in the properties of Rydberg ions of alkali earth metals [20] in 
an electric quadrupole trap. They possess a single valence electron with the remaining 
electrons forming closed shells pQ. For the description of such system one can employ 
an effective two-body approach in which the ion is modeled by a two-fold positively 
charged core (mass M, position r c ) and the valence electron (mass m, position r e ). The 
corresponding interaction (model-)potential depends only on the relative coordinate 
r c — r e and also on the angular momentum state of the atom. The latter dependence 
accounts for the quantum defect - a lowering in the energy for low angular momentum 
states in which the valence electron probes the inner electronic shells. High angular 
momentum states (typically with angular momentum quantum number I > 5) do not 
exhibit a significant quantum defect since here the valence electron is located far away 
from the ionic core thus experiencing a bare Coulomb potential. 

We now formulate the Hamiltonian of a single Rydberg ion in the presence of the 
electric potential of the Paul trap. We add a linear potential, corresponding to a time- 
dependent homogeneous electric field, $/(r) = f (t) ■ r, such that the combined electric 
potential reads $(r, t) = $r(r, t) + $;(r, t). Below, we will employ these additional MW 
fields to electronically couple different Rydberg states. This will allow us to generate 
large oscillating dipole moments, which give rise to strong dipole-dipole interactions 
among dressed Rydberg ions. 




(1) 
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Writing the interaction potential between the valence electron and the atomic core 
as V(\r e — r c |) and taking into account the coupling of the individual charges to the 
electric potentials we find 

2 2 

#iab = ^7 + 7T- + V{\r e - r c |) + 2e$(r c , t) - e$(r e , t) + H FS . (2) 
2M 2m 

The term H FS accounts for the spin-orbit coupling giving rise to the fine-structure of 

electronic levels. Its effect will be discussed in the following subsection. 

We introduce center of mass (CM) R and relative coordinates r according to 

m M 
r c = R— — r , r e = R+ — r. (3) 

M + m ' M + m v ; 

Exploiting that the nuclear mass is much larger than the electronic one and hence 
M 3> m we have r c « R and r e « R + r. Within this approximation the Hamiltonian 

becomes 

H = + + v (|r|) + 2e$(R) t] ~ e$(R + r ' t] + Hfs (4) 

P 2 p 2 



+ e 



^ 5 $(R,t) 1 <9 2 $(R,t) 

$(R ' t} 5R-- r -25^^xT^. 



FS- 



Corrections to this Hamiltonian scale with m/M which is typically about 10~ 5 . 

In case of ions in low lying electronic states the potential ([T]) provides static 
confinement along the longitudinal (z-) direction. However, transversally at no instant of 
time a confining potential is present. One rather finds a periodically oscillating potential 
saddle centered at the origin of the coordinate system. Due to the rapid periodic 
change of the confining and non-confining direction, however, the ions experience a 
ponderomotive potential that can provide transversal confinement [21] . In order to 
make this manifest we transform into a frame which oscillates at the RF frequency u in 
the CM coordinate system. This is achieved by the unitary transformation 

(pry r -i \ 

-%— \X 2 - Y 2 \ sin(u;t)J . (5) 

By applying this transformation to Hamiltonian (jSJ) one obtains 

8U 

H' = UHU^ + ih—rf = H CM + H cl + # C M_ci + H mm (6) 



with 



#cm = ^ + \mu 2 z Z 2 + \Mj p (x 2 + Y 2 ) (7) 

H el =^- + V (\r\)-e^(r,t)+H FS (8) 

Hcm-ei = -2e[a cos(wt) (Xx -Yy)-(5 (Xx + Yy- 2Zz)} (9) 
2ect 
~Muj 

— - (X 2 + Y 2 ) cos(2wt) + ef (t) R. (10) 



H m m = ~ — sin(^t) (XP X - YP y ) 
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Here Hqm provides harmonic axial and transversal confinement of the CM motion with 

the corresponding trap frequencies u z = 2-y/|f and tu p = y/2\J (jj^ — §§, which are of 
the order of a few MHz and satisfy u p ^> u z . For Ca + ions, an RF frequency uo = 2tt x 15 
MHz and the gradient parameters a = 10 9 V/m 2 and (3 = 10 7 V/m 2 the axial and radial 
trap frequencies evaluate to uj z = 2n x 1.56 MHz and u p = 2ir x 5.64 MHz, respectively. 

The term H e \ contains all dependencies on the electronic coordinates describing the 
motion of an electron in the field of a doubly charged ionic core which is superimposed 
by the electric potential $(r, t). The electronic dynamics takes place on a much faster 
time scale compared to the CM motion in the trap. The external electric field prevents, 
unlike in the field-free case, the separation of the CM and relative dynamics. The 
coupling between these motions is accounted for by Hcu-ei- Due to the large separation 
of time scales of electronic and external dynamics we will treat this coupling within an 
Born-Oppenheimer approach below. 

Finally, we have the term H mm which gives rise to the micromotion causing a 
coupling between the static oscillator levels of H CM . It can be shown [21] that for large 
enough values of the RF frequency u this coupling can be neglected and the external 
motion of the ions can be considered as if it was taking place in a static harmonic 
potential. The effect of the additional micromotion term ef (t) ■ R arising from the MW 
dressing fields can likewise be neglected in the following, since typical MW frequencies 
are of the order of at least one GHz and therefore far from being resonant with the 
external trapping frequency of the ions. 

With the full Hamiltonian ([6]) for the internal and external dynamics and the 
coupling among them at hand, we are now in the position to analyze the electronic 
properties of a trapped Rydberg ion and the mutual interplay of internal and external 
dynamics. This will be addressed in the next two subsections. 



2.2. Electronic Properties 



We proceed by inspecting the electronic Hamiltonian H e \ in Eq. (jSj), 
H el = ^ + V(\r\) + H FS + H ci 



(11) 



with the electric trapping and MW dressing fields contained in 
H e f = — e$(r, t) = H stsit + H osc + H MW 



ep x 2 + y 2 - 2z' 



2 2 

x -y 



ef(t) • r. 



(12) 



— e a cos (ut) 

In the absence of electric fields the ionic Rydberg states can be classified by the principal 
quantum number n, the angular momentum quantum number I, the total angular 
momentum j and its magnetic quantum number m. The quantum states are represented 
by \n, l,j, m) = \n,l) \j,m) a which factors in the radial part \n,l) and the angular 
momentum part \j, m) . The latter is constituted by a linear combination of products 
of the Spherical Harmonics and spin orbitals. The corresponding energies of the Rydberg 
levels are given by the well-known formula E n ij = — 4E^ y d/(n — S(l)) 2 + Eps(n,l,j) 
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Figure 2. a: Sketch of the Rydberg level structure of the Ca ion in the field free 
case without finestructure. States with I > 5 do not exhibit a quantum defect and are 
located in a degenerate manifold of states. States with low angular momentum are split 
away from this manifold. The corresponding quantum defects are provided (taken from 
Ref. [20). b: Levels of the s and p manifold in the presence of the electric potential 
<£(r, t). Its static components lead to a first order shift (quadrupole shift) of the energy 
levels according to eqs. (fT5|) and HH). In general the electronic energy is lowered 
with respect to the field- free value (sketched by the solid gray lines). The shifted 
levels are coupled by the microwave Hamiltonian Hmw (blue, solid lines) and the 
oscillating components of the trapping field (red, dashed lines). We consider the limit 
of large finestructure splitting in which the coupling between the states \n,p, 1/2, m) 
and \n,p, 3/2, m) (thin dashed lines), which is caused by H osc , can be neglected. 



where En yd = 13.6 eV is the Rydberg constant and 5(1) the quantum defect [TJ [20] . 
The energy Eps(n,l,j) accounts for the finestructure splitting due to Hps- The typical 
Rydberg level structure (without finestructure) is sketched in Fig. [2^ for the case of 
Ca+. 

In this paper we focus on states with large quantum defect, i.e. s and p states. This 
is motivated by the fact that these states are most easily accessible via laser excitations 
from the electronic ground states, and that they are energetically far separated from 
the degenerate manifold of states of higher angular momentum states. Generically, the 
energy separation AEij + i between these states scales as n~ 3 for large n. In case of Ca + 
the energy separation between the s and the p level is AE S}P /h ~ 280 GHz for n = 60. 

Let us now inspect the effect of the electron field interaction H e f contained in the 
electronic Hamiltonian (Tl2l) . We assume that the finestructure splitting of the p state 
is sufficiently large such that neither the oscillating nor the static part of H e { cause 
significant coupling between the states \n,p, 1/2, m) and \n,p, 3/2, m) and j remains 
a good quantum number. Since the finestructure splitting scales proportional to n~ 3 
and becomes larger with increasing atomic mass the validity of this assumption can be 
ensured by either choosing not too high principal quantum numbers or by employing 
heavy ions. Moreover, we can neglect the coupling between the s and p states which is 
justified by their typical energy splitting of several 100 GHz. Considering first the static 
part -f/stat °f Eq. (|T2|) we find that it gives rise to the following quadrupole shifts which 
depend on j and m: 

AE S = (13) 
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r 2 n,p) 


2 m — j 

l l J 2 



(14) 



Here (n, / 



n 



I) denotes the radial matrix element of r 2 calculated using the radial 
eigenfunctions that belong to the atomic interaction potential V(|r|). The total energy 
of the electronic states is hence given by E n ij m = E ni j + AEij\ m \. The quadrupole shifts 
with respect to the unperturbed energy of the p states are sketched in Fig. [2b. 

For too strong electric field gradients there is the danger of ionizing the ion due to 
if s t a t- To estimate the ionization threshold we consider the potential 

2e 2 
47reor 

Here we have approximated the interaction of the ionic core with the valence electron 
by a pure Coulomb potential. V possesses two saddle points located at the z-axis at 



V 



+ e/3 



2 , 2 

x + y 



2z' 



(15) 



•2sad 



± 



47reo/3 



1/3 



where it assumes the value V^ ad 



e 5 /3 



1/3 



Solving E Ry 



^sad 



yields an estimate for the gradient (3 at which field ionization would occur classically, 



7 3 

em 



1.44 x 10 21 — x 

m 2 n b 



V 



(16) 



' 27 (47re ) 4 ft 6 ri 6 

For n = 50, for example, we find (3- loll = 9.2 x 10 10 V/m 2 which is one order of magnitude 
larger than the highest gradients actually used in experiment [12] . Hence, up to n = 50 
there is no danger for the Rydberg ion to undergo field ionization. 

Let us now turn to the oscillating part of the electric field which is accounted 
for by the term H osc in the Hamiltonian ffl2|) . Since the oscillation frequency u is 
typically 2tt x 10. ..100 MHz, it is not sufficient to yield a resonant coupling between 
the s and the p states or between the fine structure components of the p manifold. 
However, by estimating | (n,p \ r 2 | n,p) \ ~ a^n 4 , where ao is Bohr's radius, we find 
that (H osc ) ps e a a^n 4 cos(u; t) and hence, although no resonant coupling is present, 
the strong modulation amplitude, which grows proportional to the fourth power of 
the principal quantum number n, might give rise to significant level shifts. From the 
symmetry properties of H osc one concludes that only levels whose magnetic quantum 
number m differs by 2 are coupled. These couplings are indicated in Fig. [2Jd by the red, 
dashed arrows. Since we assume that the finestructure splitting is much larger than 
| (Hosc) Imax the j = 1/2 states are unaffected by H osc . By contrast strong effects are to 
be expected in the j = 3/2 manifold. However, in what follows we will exclusively work 
with the Rydberg states \n,s) = \n,s, 1/2, m) and \n,p) = \n,p, 1/2, m) (see Fig. [2]). 

Finally, the term i?MW m Eq. (TT2I) . which accounts for the additional MW dressing 
fields, is to be discussed. This will be subject of the following subsection. 



2.3. Microwave Dressing of Rydberg Levels 

In this section we describe how MW dressing fields can be employed to generate strong 
interactions in an ion chain. As will be shown in Sec. 13.11 Rydberg ions aligned in a 
Wigner crystal in the Paul trap do not exhibit permanent dipole moments, and residual 
Van-der-Waals interactions are small. Thus, our aim is to generate strong interactions 
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Figure 3. Microwave dressing of ionic Rydberg states. A linearly polarized 
Dichromatic microwave field is used to couple one s with two p levels one of which 
is far detuned. After an adiabatic elimination of the latter one obtains an effective 
two-level system. This dressing is used to tailor the interaction between Rydberg ions. 



using (near-) resonant MW dressing fields to couple electronic s and p Rydberg states 
as indicated in Fig. [2b. We show in the following that the MW dressing gives rise to 
large oscillating dipole moments, which in turn lead to strong dipole-dipole interaction 
between Rydberg ions. 

The energy separation between s and p Rydberg states belonging to the same 
principal quantum number and thus as well the typical frequencies of the MW dressing 
fields are typically in the order of a few hundred GHz (cf. Sec. 12.21) . Since the 
trapping frequencies determining the external motion of the ions are in the MHz range 
(cf. Sec. I2.ip . coupling to the external motion is negligible and the MW fields exclusively 
affect the electronic degrees of freedom. Ideas similar to the MW dressing scheme 
described below were applied in the context of cold polar molecules, where a combination 
of static electric and MW fields was used in order to tune intermolecular two- and three- 
body interactions [22]. 

In our setup we apply a linearly polarized MW such that non-zero transition dipole 
matrix elements occur between the states |n, s) and \n',p). For our purposes we choose 
a bi chromatic microwave- field of the form f (t) = E\ e z cosc^it + E 2 e 2 coscu 2 t such that 
the coupling Hamiltonian reads 

Hmw — —d z [Ei coslg>i£ + E 2 cosu^t] (17) 

where d z = e z is the operator of the dipole moment. We consider three levels (one s level 
and two p levels) which are coupled by Hmw as depicted in Fig. [31 The frequencies iO\ 
and uj 2 bridge the energy separations E\ n ^ — E\ n i >p \ and -E| n , P ) — E\ H;S \, respectively, and 
are assumed to differ significantly. This can be achieved by choosing not too close values 
of n and n' . In the rotating frame and after performing the rotating wave approximation 
the electronic Hamiltonian of a single ion reads 

#3ievcis = h/\ 1 \n',p){n',p\ - HA 2 \n,p)(n,p\ 

+ - [(hQ 1 \n,s)(n',p\ + hQ 2 \n,p)(n,s\) + h.c] . (18) 

Here we have introduced the microwave detunings A x = lo\ — (-B| n ,s) — ^|n»)/^ an d 



Trapped Rydberg Ions: From Spin Chains to Fast Quantum Gates 



10 



A 2 = uj 2 — (E\ n>p \ — E\ n 3 \)/h and the Rabi frequencies f2 li2 = which involve 

the radial dipole matrix elements d\ = e(n, s\r\n',p) and d 2 = e(n, s\r\n,p). We assume 
that the MW field with frequency uj\ is far-detuned and only weakly couples the states 
\n',p) and \n,s). Furthermore it is assumed that |Ax| 3> l^l? |A 2 |- Under the condition 
Qi <C | Ai| we can adiabatically eliminate the \n',p) state and obtain an effective two- 
level system consisting of the \n,p) state and a dressed state (see also Fig. [3]) 

\ s ') = \ n ,s) ~ ^-\n',p) (19) 

with r] = fii/(2Ai), 77 <S 1. This effective two-level system can be mapped onto a spin- 
1/2 particle by identifying the states \s') and \n,p) as eigenstates of the spin operator 
5* 2 with positive (negative) eigenvalue. The Hamiltonian ffl8|) then reduces to 

with an effective magnetic field h = (Q 2 , 0, A' 2 ), detuning A' 2 = A 2 — flf/ (4A X ) and the 
spin operator S = (S x , S y , S z ). 

Due to the weak admixture of the state \n',p) the dressed state \s') obtains an 
oscillating dipole moment. The matrix representation of the dipole operator d z in the 
set of states \s'), \n,p) is given by 

a - I ( -fthcoBUit dae-*** \ 

dz ~ 3\ d 2 e^ ) ' 1 j 

Using the abbreviation D\ 2 = {di^ 2 /3) 2 one obtains the following representation of the 
dipole moment operator: 

d z = -ry/5^cos(o; 1 t)(l + 2S Z ) + 2^Jr7 2 (cos(u 2 t)S x + sm(uj 2 t)S y ). (22) 

The magnitude of the induced dipole moments is determined by the transition dipole 
matrix elements c?i j2 which can be roughly estimated as d\ )2 ~ ea^n 2 . Thus the dipole- 
dipole interaction energy scales ~ n 4 for large n. This is to be compared with the 
radiative life time of Rydberg states, which for large n and low / scales according to 
~ n 3 [TJ, thereby favoring larger values of n. We return to the question of radiative decay 



and the validity of our analysis in Sec. I3.2[ where we will use the representation ( 1221) 
of the dipole moment operator to derive an effective spin chain Hamiltonian describing 
the dynamics of Rydberg excitations in the ion chain. 



2.4- Coupling Between Internal and External Dynamics 

In this section we will analyze the effects of the coupling of internal and external degrees 
of freedom described by ifcM-ei in the Hamiltonian ([6]). We treat these coupling terms 
using a Born-Oppenheimer approach, which is reasonable since the electronic dynamics 
takes place on a much faster time scale than the external motion of the ions. In this 
framework we treat the CM coordinates as parameters while diagonalizing the electronic 
Hamiltonian. We evaluate the energy of the \n, s)-level by second order perturbation 
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theory considering only the coupling to the next |n,p)-level. After averaging over one 
RF cycle one obtains the energy correction 



e(X, Y, Z) = E M + \mJ z 2 Z 2 + \mJ 2 {X 2 + Y 2 ) (23) 



with 



32 e 2 /? 2 



J >= ~3^At; K "' s|r| "' p)|2 ' (24) 

and the energy difference AE StP = E\ n , p ) — E\ UiS ). Hence, ions excited to Rydberg 
states experience modified transversal and longitudinal trap frequencies in comparison 
to their ground state counterparts. For the s state under consideration the trap becomes 
shallower and the new trap frequencies are given by u P:Z = u p ^ z Jl + {oj' pz /uj p ^) 2 ~ 
u p ,z + uj' p 2 J (2uj p>z ) = lu PjZ + 5lu P)Z , yielding 

LO z 6 £^&s,p 

6ujp = 1 a 2 + 2p 2 \(n,s \r \n,p) r \ 2 

u p Qa 2 /{Mu 2 )- (3/e AE S)P ' 1 ' 

By estimating \(n, s \ r \ n,p)\ ~ a^r? one finds that the frequency shift 5u PtZ scales 
proportional to n 7 . For typical trap parameters, i.e. a>/J, the modification of the axial 
trap frequency is much smaller than the corresponding one for the radial frequency. For 
n = 50 and the parameters presented in Sec. 12.21 we find (5u> p )/ujp = 3.5 x 1CT 2 and 
{8uj z )/uj z = 7.4 x 10~ 4 . Below we are interested in the regime where these modifications 
of the trapping frequencies are negligible. We remark, though, that due to the scaling 
~ n 7 they can become significant even for moderately larger principle quantum numbers. 



3. Interacting Trapped Rydberg Ions 

We now turn to the discussion of the interaction between several ions. First, we 
analyze the various types of emerging interactions between Rydberg ions, and discuss 
the interplay of interionic interactions and the effect of the external trapping fields. 
Then we proceed with the discussion of interacting Rydberg ions dressed by MW fields 
and derive the corresponding effective spin- 1/2 Hamiltonian describing the Rydberg 
excitation dynamics in the ion chain. The role of the two effective spin states of each 
ion will be played by two Rydberg states, as outlined in Sec. I2.3l and illustrated in Fig. [3j 
We illustrate the effective spin dynamics by studying the process of resonant 
excitation transfer in a mesoscopic chain of ten ions. Furthermore, we suggest a 
scheme, which allows to observe the effective spin dynamics in dressed ground state 
ions. We finally show how the trapped Rydberg ion system can be exploited in the 
context of quantum information processing and discuss an implementation of a two- 
qubit conditional phase gate. 
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Figure 4. Interacting Rydberg ions. The net charge of the ions leads, apart from the 
common Coulomb repulsion, to a charge-dipole and a charge-quadrupole interaction 
both of which are absent in systems of neutral Rydberg atoms. 



3.1. Static Interactions between Trapped Rydberg Ions 

We consider the interaction between ions % and j with coordinates given by (Rj, r^) and 

(Rj,rj) as depicted in Fig. [U The Coulomb interaction Vfa, rj, Rj, Rj) between the 

charges belonging to different ions can be written as 

VCRi^Tj,^) 4 2 

e 2 /(47re ) |R>-R;| |Ri - (R,- + r 3 -)| 

2 1 

~ |(R. ; + ri) - Rj\ + |(Ri + n) - (Rj + rj)!' (28) 
We assume that |Rj — Rj| ^> which is very well fulfilled, since - as discussed above 
- we are interested in the parameter regime where the average interparticle distance ( 
in the ion trap is much larger than the extension of the electronic wave function OR y of 
the Rydberg ions. Performing a multipole expansion up to second order in the small 
parameter (oR y /C) an d abbreviating |Rj — Rj| = |Ry| = Rij and = (Rj — Hj)/Rij 
we obtain the following form for the Rydberg ion-Rydberg ion interaction potential: 
V(Ri, Rj, rj, _ 1 (Rj - Rj)(r< - r,-) 



e 2 /(47r eo ) Rij ' R% 



, r 2 -3(n tJ T t ) 2 + r|-3(n tJ - rj ) 2 
H 



1R% 



rj • - 3(rij^rj)(nj i • r jy 
R 3 

The first term accounts for the Coulomb interaction between two singly charged ions 
and is independent of the degree of electronic excitation. In case of Rydberg ions the 
displacement of the electronic charge from the ionic core leads to the exhibition of a 
dipole moment which interacts with the charge of the other ion. This dipole-charge 
interaction gives rise to the second term. The third term accounts for the charge- 
quadrupole interaction. These three terms are absent in the case of interacting neutral 
Rydberg atoms. The last term is the well-known dipole-dipole interaction, which is also 
present in neutral systems. 

For iV ions stored in a Paul trap, at sufficiently low temperature and tight radial 
trapping, u p ^> lu z , the ions form a one-dimensional Wigner crystal with equilibrium 
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positions along the Z-axis [TJ)J [T5], which are determined by the interplay between the 
Coulomb repulsion among the ions and the external confinement by the trapping fields. 



In Appendix A we show that after a harmonic expansion of the Hamiltonian around the 
equilibrium positions zf^ (given by Eq. (1A.3|) ) the full Hamiltonian of A" interacting 
ions can be written as 

N 

Hi on s = i?ph + J! ^el,i + -^int-ext + -^dd- (30) 

i 

The first term 

N 

-^ph = ^a,n%L,n a <x,n (31) 

a=x,y,z n 

describes the external oscillation dynamics of the ionic cores with a^ a n and a a ^ n being 
the respective creation and annihilation operators of the normal modes (cf. Eqs. (1A.5I) - 
( 1A.8j) ). The second term determines the electronic level structure of the ions: the 
charge-quadrupole term gives rise to a position dependent variation of the electric field 
experienced by the trapped ions, which can be absorbed in the single particle ion-field 
interaction H e f given by Eq. (TP2]) . Thus, the electronic Hamiltonian of the i-th ion takes 
the form 



2 



#ei,i = p- + V{\T i \) + e/3' i xl + yj- 2z] - e a cos (cut) 
2m 1 J 



2 2 



-ef(t)-Ti (32) 
with position dependent gradient 

e N 1 

tt = e+^- o E lzm _ zm =/»+w- (33) 

The third term ifi nt _ ex t in Eq. fl30|) accounts for the coupling of the internal and external 
dynamics, which is partly due to the dipole-charge interaction and partly due to the 
inhomogeneous electric field of the Paul trap (see Eq. (1A.12I) ). The resulting coupling 
between the internal and external dynamics becomes also position dependent and thus 
leads to a state-dependent variation of the trapping frequency (see Sec. 12.41) . However, 
since we work in a regime where these shifts are negligible already in the single-ion case 
they can also be neglected in the case of many ions. 

Finally we have the dipole-dipole interaction in Eq. (130!) . which after the harmonic 
expansion is given by 

The discussion of the individual terms contained in Hamiltonian (130!) proceeds in 
analogy to the one presented for the single ion case. The electronic dynamics of each ion 
is governed by H c \ i in Eq. fl30l and is therefore dependent on the equilibrium positions 
of the ions. This gives rise to position dependent electronic energies, i.e. E n ij m —> E n ij mti 
according to Eqs. (fT2l - (fl4l) with the ion-dependent gradients of Eq. (J33]). 
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We remark that even ions located at the edges of the ion chain experience the 
same electric field as an ion in the center of the trap. This is due to the fact that 
the electric field of the Paul trap is compensated by the field created by the ions 
themselves. Therefore, the ions do not exhibit permanent dipole moments and each 
ion can be considered to be located in a center of a local electric quadrupole field. This 
is the reason why the electronic properties of the system of trapped ions are virtually 
unaffected by the actual number of particles at hand. 

The term of interest to create strong interactions among the ions is the dipole-dipole 
interaction. In the parameter regime where the interparticle distance ( is much larger 
than the extension of the electronic wave function aR y of a Rydberg ion, the dipole-dipole 
interaction can be treated perturbatively. For the considered s and p states the expected 
interaction energy can be estimated by -EVdw ~ e 4 | (n, s \ r \ n,p) | 4 /((47reo) 2 | AE StP \( 6 ). 
For Ca + , n = 50 and ( = 5 //m one finds -E v dw ~ h x 200 kHz. In the next subsection 
we demonstrate that the MW fields introduced in Sec. 12.31 in order to manipulate the 
electronic level structure can give rise to much stronger interactions among the ions. 



3.2. Interaction Between MW-Dressed Rydberg Ions 

To study the interaction between MW-dressed Rydberg ions we first transform the N- 
ion Hamiltonian ( 130]) to a frame of reference, which rotates at the microwave frequencies 
Ui and u>2 of the dressing fields (in analogy to the single-ion transformation in Sec. 12. 3j) . 
In the rotating frame, the terms in H- mt - ex t oscillate rapidly at the frequencies uj\ and 
uj 2 and u)\ — uj 2 . Hence, in the rotating- wave approximation the internal and external 
dynamics decouple (cf. discussion in Sec. 12. ip . 

We use the dressed ion Hamiltonian (120]) and the expression (l2Tj) for the electronic 
dipole moment in order to represent the Hamiltonian H mt = Yli H e \^ + H^ which is part 
of the full Hamiltonian (1301) . The terms H e \^ are substituted by the two-dimensional 
representation given by Eq. (1201) . In this representation we can write the dipole-dipole 
interaction as 

1 N dWd(?) 

Hdd = ~^o^\Z^-Z^ (35) 



j 



where d$ is given by Eq. (T2T1) with the index i labeling the respective ion. In the 
rotating wave approximation we obtain 

N N rl -, 

H M ='Eb®8® + D l £ jvf + ViVjS® (36) 

N 

+ £ Vii [D 2 (S2>S<P + S®Sjp) + TwDiSWS®] , 

with Vij = -2/(47re |^ (0) - zf\ 3 ) and an effective magnetic field h« = h(tt 2 ,0, A ( 2 i] ). 
The coefficients Di >2 and rji characterize the coupling strengths and depend on the 
transition matrix elements between the involved Rydberg states and the MW dressing. 
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Figure 5. a: Schematic level scheme of the three MW-coupled Rydberg levels, for a 
chain of 5 ions. The individual, position-dependent energy shifts of the I = 1 states 
give rise to an inhomogeneous distribution of MW detunings (not true to scale), b: 
Schematic spin dynamics of a chain with the initial configuration: first ion in the \n,p) 
state, all other ions in the \s') state (blue filled circles). After a certain time the 
Rydberg excitation has travelled to the right end of the chain (red open circles). A 
numerical example for this excitation transfer is shown in Fig. [6] 



As depicted in Fig. [3] the effective spin corresponds to a two- level system constituted by 
two MW-dressed Rydberg states (see Sec. 12.31 for details). 

The terms linear in the spin operators represent a coupling of a series of spins to 
an inhomogeneous effective magnetic field whose strength and direction are determined 
by the position dependent electronic energies of the Hamiltonian (!32l) . The terms in 
the second line of Eq. (1361) (quadratic in the spin operators) represent a ferromagnetic 
Heisenberg chain with 1/r 3 exchange-type interaction [23]. Establishing the connection 
to neutral Rydberg gases the S^S^ and S^S^ can be interpreted as resonant dipole- 
dipole interaction terms. The term which is proportional to S^S^ resembles the 
interaction of two static dipoles. In general, the coefficients rji = f2^/(2A^) depend on 
the ion index, since different energy shifts of the \n',p) level due to position-dependent 
charge-quadrupole terms (cf. Eqs. (|32|) and (|33|) ) lead to different detunings for 
different ions (see Fig. [5]). In case one does not admix the \n',p) state to \n, s) (i.e. Qi = 
in Eq. (|T8|) ). the coefficients rji vanish and the interaction Hamiltonian solely describes 
resonant dipole-dipole interaction in the presence of an effective magnetic field with a 
constant component in x and a position-dependent component in z-direction. 

The physical realization of effective spin dynamics, as provided by the Hamiltonian 
( 1361) , has been of significant interest in atomic physics during the last few years as "analog 
quantum simulators" of (mesoscopic) condensed matter systems. The distinguishing 
feature of the present setup is the large coupling strength between effective spins, which 
scales proportional to n A and is of the order h x 500 MHz (n « 50, typical interparticle 
spacing £ « 5 /mi). The characteristic time scale of the effective spin dynamics is thus of 
the order of a few nanoseconds. This is significantly shorter than the typical decoherence 
time in the system, which is set by the radiative lifetimes of the involved Rydberg states, 
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which scale as n 3 and are typically of the order of /is (for Ca + ions and n = 50 the 
lifetime is ~ 10 /xs [20]). 

Realization of effective spin models has also been proposed in the context of trapped 
ions in their electronic ground state, and with cold atoms and polar molecules in 
optical lattices. For trapped ions involving electronic ground states models analogous 
to fl36l) can be derived where the typical coupling strengths for the effective spin-spin 
interactions are in the range of tens of kHz [16] with (long) decoherence times as 
described in the context of ion trap quantum computing [2U EH]. Effective Heisenberg 
models with nearest neighbor interactions are also obtained with cold atoms in optical 
lattices [26, 27J, where the time scales of exchange interactions can be of the order of a 
few hundred Hz [28[ |29]. We note that these energy scales are also directly related to the 
temperature requirements for the preparation of an effective zero temperature ensemble. 
Finally, effective spin models, such as the Kitaev model |30j, have been proposed with 
polar molecules in optical lattices [31j [321 [33] . In this context electric dipole moments of 
a few Debye can be induced by external DC and microwave electric fields in the rotational 
manifold of molecules prepared in their electronic and vibrational ground state. The 
electric dipole-dipole interactions, which are strong and long range in comparison with 
neutral atom collisional interactions, can lead to effective offsite-coupling strengths of 
the order of tens or hundred kHz, limited mainly by the conditions imposed by optical 
trapping, while coherence times are of the order of seconds, as determined e.g. by 
spontaneous emission in off-resonant light fields forming the trapping potential. 

A second feature of simulating spin dynamics according to Hamiltonians of the type 
( 136]) is the single ion addressability and read out [M] which the present setup inherits 
from the experimental developments in ion trap quantum computing. In contrast, 
neutral atoms and molecules in optical molecules ususally allow global addressing by 
laser light, even though a significant effort is devoted at the moment to develop these 
tools also for optical lattice setups j35[ [36]. Note, however, that neutral atoms and 
molecules in optical lattices will typically allow for systems with a significantly larger 
number of "spins" than in the ion case. 

One of the experimentally most challenging aspects of realizing a Rydberg ion chain 
is the requirement of 7r-pulses to transfer ions from the ground state to the Rydberg 
state. In Sec. 13.41 we discuss a version of an effective spin chain where the Rydberg 
dipoles are admixed to the electronic ground states with an off-resonant laser process, 
resulting in ground state ions with effective oscillating dipole moments. 

3.3. Resonant Excitation Transfer 

As an illustration of the spin dynamics contained in Hamiltonian (136]) . we study the 
transfer of an excitation from one side of the Rydberg ion chain to the other end. For 
simplicity we choose a scenario in which rji = (Q\ = in Eq. (fl8|) ). As discussed in 
the previous subsection in this case the interaction reduces to purely resonant dipole- 
dipole interaction. The ion-dependent coefficient of the dipole-dipole interaction can be 
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Figure 6. Transport of an excitation through a chain of ten ions. The ion at site 
1 is initially excited to the state ||) while all others are in the state |!). After the 
time t = 1.8 Ti/ J the excitation is transferred from the first to the tenth ion. The inset 
shows the time evolution of the expectation values (Sz ) and (Si 10 ^). 



written as 



2Mut 



47re n Z 



(0) 



'(°)|3 



\Uj — Uj 



(37) 



where the Ui are the equilibrium positions in dimensionless coordinates, Ui = Z^/( 

-I/O 

with ( = [e 2 1 (4ire Mu! 2 )} (see Ref. [15J). The scale of the interaction energy is then 
given by J = — 2M e 2 z D 2 . We choose the MW frequency u 2 such that it is on resonance 
with the energy gap between the levels \n, s) and \n, p) determined by the gradient ( 1331 ) 
with 5 /3i = 0. Following Eq. (fT4"|) . the position dependent change of the gradient 5 fa 



gives rise to a position dependent variation of the detuning which reads 

4 



A 



2,i 



edfa (n,p \ r 2 \ n,p^} 



n,p \ r | n,p 



N 

£ 

3^ 



Ui 



Ui 



N 



Ui 



3 ' 



(38) 



This situation is depicted in Fig. [U The effective magnetic field B z and the coupling 
constant J do not scale independently since both D 2 /e 2 and the matrix element 
(n,p | r 2 | n,p) are of the order of a\n^. The precise value depends on the ionic species. 
For our simulations we use the parameters B z = 0.65 J and hQ 2 = 0.01 J. Initially 
the system is prepared such that the first ion (k = 1) is in the state ||) = \n, s) while 
all others are in the state ||) = \n,p). This state is sketched in Fig. [5b by the solid 
circles. The temporal evolution under the Hamiltonian fl36l) leads to a transfer of the 
excitation from the first ion to the last ion in the chain (open circles in Fig. [5b)- The 
corresponding numerical data is shown in Fig. [6] where we monitor the time evolution 
of the expectation value (S]p) in a chain of 10 ions. The excitation transfer from the 
first to the tenth ion takes place in a time t = 1.8H/J with an efficiency of 89 %. Since 
J can be in the order of H x 500 MHz this resonant excitation transfer can therefore be 
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Figure 7. MW dressing of ground state ions. Two lasers with Rabi frequencies f2 s , 
rip and detunings A s , A p weakly couple two electronic ground states \gi), \g2) to 
the Rydberg levels \n,s) and \n,p), respectively Thus ions residing in the electronic 
ground states are dressed and obtain time-dependent oscillating dipole moments. 



achieved in about 3.6 ns which is much less than the lifetime of ionic Rydberg states. 
3.4- Spin Dynamics in the Electronic Ground States 

In view of the short transition wave length of 100...125nm associated to a transition 
from an electronic ground state to a Rydberg level [37] , it is experimentally challenging 
to realize 7r laser pulses, which transfer the entire electronic population to the Rydberg 
levels, as required for the initialization step of the chain of effective spins. Thus, we 
outline an alternative scheme, which does not require the transfer of the full electronic 
population to Rydberg states and which is based on an adiabatic admixture of Rydberg 
levels to two electronic ground states by near-resonant CW lasers. 

We extend the 3-level scheme of Sec. 12.31 by including two ground states \g\) and 
|<72), which are coupled by two lasers to the (undressed) Rydberg states \n, s) and \n,p), 
respectively (see Fig. [7J. As before, we use the MW dressing fields to couple the set of 
bare Rydberg states \n',p), \n, s) and \n,p). The two CW lasers weakly couple the two 
ground states to the resulting dressed Rydberg level structure. Weak coupling requires 
that the laser Rabi frequencies are small compared to the MW Rabi frequencies and 
that the laser fields are sufficiently far detuned from the three dressed Rydberg states. 
Due to the laser coupling the two dressed ground states \g[), \g' 2 ) obtain time-dependent 
oscillating dipole moments, which finally lead to dipole-dipole interactions between laser- 
dressed ground state ions. The dressed ground states now constitute the two-level system 
of interest (effective spin degree of freedom) and play the role of the dressed Rydberg 
levels \s') and \n,p) of Sec. 12.31 



As shown in |Appendix B| for a specific set of Rabi frequencies and detunings of the 
MW and laser fields the representation of the dipole operator in this effective two-level 
system reads 



gL - 1 



c v .\c s ^\fFh.^ 2t 2c Si2 c p ' i 2V / A"cos(^it) 
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(39) 



with the factors c a given by Eqs. (IB. 71) . Thus, the resulting dipole operator as well 
as the resulting spin chain Hamiltonian describing the dynamics in the dressed ground 
states \g[), \g' 2 ) have the same structure as in Eqs. (1221) and ( 1361) . with the difference that 
the magnitude of the dipole moments is decreased by the factors c P) iC Sj2 and c p ' j2 c S)2 , 
respectively. 

Due to the weak admixture of the Rydberg states the dressed ground states 
have a finite lifetime, i.e. the time until a photon of the dressing lasers is scattered. 
Denoting the lifetime of the involved Rydberg states by r = 1/T, this scattering rate 
is approximately given by r scat = |c a | 2 r, where c a represents a typical value of the 
admixture coefficients in Eq. ( IB. 71) . For instance, the lifetime of the state \g[) is enhanced 
by a factor of |c Pi i|~ 2 with respect to the radiative lifetime of the Rydberg state \n,p), 
since only the fraction |c Pi i| 2 of the electronic population resides in the Rydberg state. 
Moreover, the fraction of the admixed Rydberg states affects the interaction among the 
ions. Compared to the bare dipole-dipole interaction the interaction between the dressed 
ions is suppressed by a factor ~ |cq,| 4 resulting in an effectively slower spin dynamics. 
The typical time scale for the spin dynamics, e.g. the excitation transfer discussed in 
Sec. 13.21 thus increases proportional to |c a |~ 4 . Hence, comparing the scaling of the 
interaction strength and the lifetime implies that the adiabatic admixture of Rydberg 
states must not be too weak in order to avoid that decoherence due to spontaneous 
emission becomes an issue during the temporal evolution of the spin dynamics. 

3.5. Fast Two-Qubit Quantum Gates with Rydberg Ions 

In most two-qubit gate schemes, which have been suggested and implemented so far with 
trapped ions (see. Refs. [381 E3 HQ] for a few examples), the motional sidebands of the 
ions have to be spectroscopically resolved, which limits the achievable gate operation 
times, typically to the order of the inverse external trapping frequency. In order to 
overcome this limitation, it has been proposed to use specifically shaped trains of off- 
resonant laser pulses to implement geometric two-qubit ion gates with much faster gate 
times [41] . In the context of quantum information processing with neutral atoms, 
potentially fast two-qubit gates can be achieved for pairs of Rydberg atoms in an 
optical lattice, based on the strong and long-range dipolar interactions among the atoms 



For the system of trapped Rydberg ions we have shown that despite the absence 
of permanent dipole moments MW-dressing fields can be used to generate strong 
dipole-dipole interactions between Rydberg ions. We now aim to exploit the electronic 
interaction Hamiltonian (|36|) for the implementation of a fast conditional two-qubit 
phase gate along the lines of the proposals developed for neutral Rydberg atoms. 
Such gate is characterized by the truth table \g a )m\9b)n e^ a_2 ^ fe_2 ^ cnt |<7 ) m |<7&) n with 
a, b = 1, 2 labeling the two ground states, and the ion indices m, n. Thus, the two-qubit 
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Figure 8. Level scheme for the implementation of a conditional two-qubit phase 
gate. The ground states \gi) of the respective ions are coupled via a laser with time- 
dependent detuning A s (t) and Rabi frequency f2 s (i) to the Rydberg level \n, s). A 
MW field of constant Rabi frequency H2 and ion-dependent detuning A2. m couples 
the Rydberg levels \n,s) and \n,p). 



state \gi)m\gi) n obtains a phase shift, while the other three states are unaffected (up to 
trivial single qubit phases) [43J. 

We identify the two ground states |<7i) m and \g2)m of each ion m as logical qubit 
states |0) and |1). The ground state \gi) is coupled to the Rydberg state \n,s) by a 
near-resonant laser with time-dependent effective Rabi frequency Q s (t) and detuning 
A s (t) = u> s (t) — (E\ n>s \ — E\ gi \)/h, which can be chosen equal for both ions m and n. 
The second ground state \g 2 ) is not coupled to any Rydberg state. We again consider 
the scenario in which only one additional MW field with Rabi frequency ^2 arid ion- 
dependent detunings A 2jm (cf. Sec. 13.21) is applied (see Fig. i.e. Oi = and rj m = in 
Eq. ( 1361) . To perform the gate operation, we apply laser pulses to the two ions (similarly 
to the adiabatic gate scheme presented in Ref. [8]). The variation of the laser pulses is 
assumed to be slow on the time-scale set by Q s and A s such that the system follows 
adiabatically the dressed states, which arise from slowly switching on the laser coupling. 
This adiabaticity condition guarantees that after applying the laser pulses the electronic 
population still completely resides in the initial electronic ground states. During the 
application of the laser pulses and the resulting dressing of the electronic ground states 
part of the electronic population is transferred from the states \gi) m and \gi) n to the 
Rydberg states, where the ions interact via the resonant dipole-dipole interaction and 
thereby accumulate an entanglement phase, given by 

0ent(t) = / dr(e| Pl )| Pl )(r) - e ]gi)]g2) (T) - e\ g2)]gi) (T)) (40) 

J 

where c\g a ) m \ gb ) n denotes the eigenenergy of the instantaneous eigenstate connected to 
the state \g a )m\9b)n in the absence of laser pulses. Fig. [9^ shows a specific choice for the 
pulse profile, the time-dependent energies of the dressed states adiabatically connected 
to the different ground states. The resulting accumulated phase shifts are presented in 
Fig. [9b. The gate operation time T is approximately two orders of magnitude larger than 
the inverse of the maximum Rabi frequency VLf^ of the applied lasers in order to satisfy 
the adiabaticity condition (for the chosen set of parameters T = 100/fi^ - 1 = 10/is). In 
order to minimize imperfections in the gate operation due to spontaneous scattering of 
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Figure 9. Two-qubit conditional phase gate between the first and second ion in a 
chain of ten ions, a: Pulse profile: Rabi frequency f2 s (t) and detuning A s (t) of the 
dressing lasers. The MW Rabi frequency is chosen fl 2 = 57.5 MHz, the ion-dependent 
MW detunings are A 2 ,i = -279 MHz, A 2 , 2 = -667 MHz. Dipole-dipole interaction 
energy scale J/h — 500 MHz (cf. Sec. 13 . 2[) . b: Accumulated phase shifts in the dressed 
states adiabatically connected to the initial ground states, and resulting entanglement 
phase 4> on t(t)- 



photons of the dressing lasers (cf. discussion in Sec. I3.4p . the gate operation has to take 
place on a time-scale much faster than the lifetime r of the involved Rydberg levels, i.e. 
T ^ r. For sufficiently high laser Rabi frequencies and long-lived Rydberg levels this 
condition can be satisfied. 

4. Conclusions and Outlook 

In this paper we have shown that trapping of Rydberg ions in a linear electric ion trap 
under realistic conditions is feasible. We have found that for not too large principal 
quantum numbers field ionization of the Rydberg ions due to the trapping fields is 
negligible, and that the coupling of electronic and external dynamics of the ions results in 
renormalized trapping frequencies for ions excited to Rydberg states. We have suggested 
to use MW dressing fields in order to generate strong dipolar interactions among the 
ions. The Rydberg excitation dynamics of the MW-dressed ions can be described by an 
effective interacting spin-1/2 Hamiltonian. The strong interactions give rise to a typical 
corresponding time scale of the order a few ns, which is substantially shorter than the 
typical decoherence time set by the radiative decay of Rydberg states. We have studied 
the dynamical transfer of a Rydberg excitation through the ion chain and discussed the 
implementation of fast two-qubit gates. While the laser excitation of ions to Rydberg 
states is an experimentally challenging task, the system offers the prospect of studying 
coherent many-body quantum dynamics on fast time scales in a well-controlled and 
structured environment. 

Beyond the present work, trapped Rydberg ions offer a rich playground for further 
studies of more involved Rydberg dynamics: Combining e.g. ions of different species 
or exciting only a certain number of the trapped ions to Rydberg states offers the 
possibility to introduce further inhomogeneities in the spatial distribution of effective 
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spins and to study in a well-controlled way spin chain dynamics in the presence of 
disorder and effective spin vacancies. The analysis can also be readily extended to study 
Rydberg excitation dynamics in two- or three-dimensional ion crystals. The invidual 
addressability of trapped Rydberg ions allows for the setup of an excitation trap [9], 
where an effective spin excitation propagating along the ion chain is trapped or extracted 
from the system once it reaches a dedicated site. 

The present work focuses on the parameter regime, where the interparticle distance 
in the ion chain is much larger than the extension of the electronic wave function of ions 
excited to Rydberg states. Another interesting scenario is the situation where those 
two length scales become comparable. If this regime can be achieved, electrons could 
hop between overlapping Rydberg orbitals of adjacent ions, which would pave the way 
towards the realization of mesoscopic Hubbard-type models in an ion trap. 
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Appendix A. Hamiltonian of N Interacting Rydberg Ions 

In this appendix we derive the Hamiltonian for a system of N interacting trapped 
Rydberg ions. Using the single ion Hamiltonian ([6]) and introducing a label for the 
respective ion we find 

N i N 

H N = Y,H[ + - £ VCH^R,-,^-). (A.l) 



with V given by Eq. ( 1281) . The external dynamics is governed by the interplay of the 
harmonic confinement and the Coulomb force between the ions. The corresponding 
potential reads 

i N i „2 N i 

V ' M = 2 M E N Z ? + "* W + >'. 2 >] + ~ E W : (A.2) 

i u 1,3 (&) V 

Following Refs. [H3, [TTj we perform a harmonic expansion of V ext which puts us 
into position to describe the external dynamics in terms of uncoupled harmonic 
oscillators/phonon modes. To this end we first calculate the ionic equilibrium positions 
Rj- " 1 . In a linear Paul trap with tight transversal confining, u p 3> uj z , the ions align 
along the z-axis of the trap such that R J -°' ) = (0, 0, zf^) where the zf^ are determined 
by (dVext/dZi) = which leads to the system of equations 



9 N 7 (0) _ 7 (0) 
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With the zf^ being determined, the harmonic approximation of the potential V ex t reads 



1 



V ext = -M J2 E K ^TQj (A.4) 



with the displacements qf = Xi, q\ = Yi and qf = Z\ — Zf' and the coefficient matrix 

f ^ ~ C «4^M ^WTiWp for i = J 

{ Ca 4ne M | Z (0)_ Z (0)| 3 ^ OT * ^ 

where c x , y = 1, c 2 = —2. The vibrational dynamics of the ID-chain is then described 
by 

i N i N 

H^ = ~M £ £K-g«g« + _ 2 E(^) 2 - (A-5) 

* «=x,j/,^ i j ^iVi a =x,y,z i 

The Hamiltonian i7 p h can be brought into diagonal form by introducing phonon modes 
via the orthogonal matrices M a 



N M a 

It = E 1 = J <n + ««,»), (A.6) 

l2Muj^ n /h 



N yra. 

n ^2/(hMu a , n ) 

N 

^M^ m = W X (A.8) 

which leads to Eq. (1311) for the phonon dynamics. 

We now proceed by expanding the charge-dipole, the charge-quadrupole and the 
dipole-dipole interaction around the equilibrium positions of the ions. The charge- 
quadruple and the dipole-dipole interaction can be approximated by 

1 p 2 N 2_o( \2 i 2 N 2 _o 2 

Cq 24vre ^ 2R% ~ 2 4vr eo 2\zf ) - zf |s 1 " J 

1 e Tj ■ Fj 3(n.jj • rj)(n.jj ■ r^) 1 e Fj • Fj 3z%Zj 
dd "24^^ i?g "24^^| Z f)_ Z (o)|3- 

(A.10) 

The charge-quadrupole term leads a position dependent variation of the electric field 
and can be absorbed in the single particle ion-field interaction if e f which is given by 
Eq. (I12p . The electronic Hamiltonian of the z-th ion then assumes the form of Eq. (1321) 
with the ion-dependent gradient ( 1331) . 

In order to treat the charge-dipole coupling we introduce the compact notation 
rf = Xi, Rf = Xi, . . . which enables us to write 

1 e 2 " (R,-R J )(r,-r,) 
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1 J2 n 7 (o) 7 (o) 

i e ^ — Zjj ^ 

^2^^\ Z ^-zfY Zl Zj) 
N e 2 1 

It turns out that by making use of the equilibrium condition (IA.3I) one can combine this 
expression with the term J2i -^CM-ei,i that arises from the summation of the single ion 
terms of the Hamiltonian ($5$). The combination then includes all couplings between the 
internal and external dynamics. It reads 

N ^ N 

# int _ cxt = -2ea cosM) x { - qV yi ] + -M £ E K ^r^ . (A.12) 

i ^ a=x,y,z i : j 

Finally, neglecting the micro-motion, the complete Hamiltonian of N Rydberg ions in 
the linear Paul trap is given by Eq. (13TT) . 



Appendix B. Effective Dipole Moment of Laser-Dressed Ground State Ions 

In this appendix we derive the form (1391) of the dipole operator for laser-dressed ground 
state ions. The Hamiltonian for the system of five coupled electronic levels as depicted 
in Fig. ([7]) in the rotating frame defined by the transformation 

U(t) = e- i »* t \g 1 )(g 1 \+e i ^-^ t \g 2 )(g 2 \ 

+ e~ iuilt \n',p)(n',p\ + \n,s)(n,s\ + e iuJ2t \n,p)(n,p\ (B.l) 

and in rotating- wave-approximation is given by H 5leve i s = H + H pert , 

H = hA s \g 1 )(g 1 \ + hA p \g 2 )(g2\ + hAx\n',p)(n',p\ - KA 2 \n,p)(n,p\ 
h 

+ - (tti\n r ,p)(n, s\ + Vt 2 \n, s)(n,p\ + h.c.) , (B.2) 

H pcvt = | (^ s \gi){n,s\ + Q p \g 2 )(n,p\ +h.c.) . (B.3) 

The MW detunings A 12 and Rabi frequencies Qi >2 are defined as in Sec. I2.3[ and the 
detunings and Rabi frequencies characterizing the two additional laser fields are given 
by A s = u s - {E ]n>s) - E lgi) )/h, A p = u p - (E ln , p) - E\ g ^)/h- A 2 and f2 s , Q p . We 
assume that the near-resonant MW fields are much stronger than the two lasers such 
that they determine the level structure of the dressed Rydberg states. Therefore we 
treat the coupling to the ground states as a perturbation H pert of the system. We now 
apply a canonical transformation to the Hamiltonian [H] , 

e~ 5 -^levels e 5 = H + H pcrt + [Hq, S] + [-ffpert, £] + ... (B-4) 

Choosing S such that [H , S] = —H peTt guarantees that to first order in the perturbation 
ifp ert the transformed Hamiltonian is block-diagonal and that the two ground states 
become decoupled from the Rydberg states. The transformation yields two dressed 
ground states, \g[ )2 ) ~ (1 + S)\gi, 2 ), 

Wi) = \9i) + c P ',i\ n ',p) + c s ,iK s) + c pA \n,p) (B.5) 
l^} = life) + c p > >2 \n',p) + c S)2 |n, s) + Cp, 2 \n,p) (B.6) 
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with 



/ 


Cp',1 N 




( -(A S + A 2 )Q 1 


\ 




C s ,l 


= 7s 


2(Ai- A S )(A S + A 2 ) 




\ 


Cp,l / 




^ {A 1 - A s )tl 2 


/ 


/ 


Cp>,2 ^ 




( -^i0 2 ^ 






C s ,2 


= 7 P 


2(Ai - A P )Q 2 




V 


C p ,2 J 




{ 4(A X - A P )A P + y 





where 



7. 



s.p 



s.p 



(A SiP + A 2 )[4(A 1 - A 
For the special choice of laser detunings 

1 



PJ s,p 



^l] + ( A s,p 



AOOI' 



we find 



-A, 



C s ,l 



V 



An 






n 2 / 



A x - JA( + Q( 



ill f2r> 



Cp',2 
Cp,2 / 



Ai+ v /a|+c| 



V 



o 



/ 



(B.7) 



(B.8) 



(B.9) 



(B.10) 



Thus, for this set of parameters, the Rydberg state \n,p) is exclusively admixed to the 
ground state \gi), while the second ground state obtains a small fraction of the Ryd- 
berg states \n',p) and \n, s). This implies that the dressed ground state \g' 2 ) possesses a 
non- vanishing dipole moment, which oscillates at the MW frequency u>i, and that there 
exists also a transition dipole matrix element between the two dressed ground states. 
We now identify the two dressed ground states \g[), \g' 2 ) with the eigenstates of the S z 
spin operator along the lines of Sec. 12.31 and obtain the dipole operator 
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